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Abstract 

It has recently been shown that triply periodic two-phase bicontinuous composites with interfaces 
that are the Schwartz primitive (P) and diamond (D) minimal surfaces are not only geometrically 
extremal but extremal for simultaneous transport of heat and electricity. The multifunctionality 
of such two-phase systems has been further established by demonstrating that they are also ex- 
tremal when a competition is set up between the effective bulk modulus and electrical (or thermal) 
conductivity of the bicontinuous composite. Here we compute the fluid permeabilities of these and 
other triply periodic bicontinuous structures at a porosity 0=1/2 using the immersed boundary 
finite volume method. The other triply periodic porous media that we study include the Schoen 
gyroid (G) minimal surface, two different pore-channel models, and an array of spherical obstacles 
arranged on the sites of a simple cubic lattice. We find that the Schwartz P porous medium has 
the largest fluid permeability among all of the six triply periodic porous media considered in this 
paper. The fluid permeabilities are shown to be inversely proportional to the corresponding specific 
surfaces for these structures. This leads to the conjecture that the maximal fluid permeability for 
a triply periodic porous medium with a simply connected pore space at a porosity = 1/2 is 
achieved by the structure that globally minimizes the specific surface. 
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I. INTRODUCTION 



Triply periodic minimal surfaces are objects of great interest to physical scientists, biol- 
ogists, and mathematicians. A minimal surface, such as a soap film, is one that is locally 
area-minimizing. Minimal surfaces are defined as surfaces with zero mean curvature. A 
remarkable class of minimal surfaces are those that are triply periodic (i.e., periodic in three 
directions). A triply periodic minimal surface is infinitely extending, has one of the crystal- 
lographic space groups as its symmetry group and, if it has no self-intersections, it partitions 
space into two disjoint but intertwining regions that are simultaneously continuous If 
one thinks of each of these disjoint regions as the phase of a composite, then those two-phase 
composites in which the interface is a triply periodic minimal surface are a special class of 
bicontinuous two-phase composites. The Schwartz primitive (P), the Schwartz diamond (D) 
and the Schoen gyroid (G) minimal surfaces partition space into two disjoint but congruent 
regions, i.e., the volume fraction of each phase is 1/2 (see Fig. Q). Triply periodic minimal 
surfaces arise in a variety of systems, including block copolymers [3], nanocomppsites j^, 
micellar materials ^, and lipid- water systems and certain cell membranes jslO, 





FIG. 1: (Color online) Unit cells of three different minimal surfaces. Left panel: Schwartz primitive 
(P) surface. Center panel: Schwartz diamond (D) surface. Right panel: Schoen gyroid (G) surface. 

One of us has very recently shown that triply periodic two-phase bicontinuous composites 
with interfaces that are the Schwartz P and D minimal surfaces are not onl y g eometrically 
extremal but extremal for simultaneous transport of heat and electricity [isi ll^ . More 
precisely, these are the optimal structures when a weighted sum of the effective thermal 
and electrical conductivities is maximized for the case in which phase 1 is a good thermal 
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conductor but poor electrical conductor and phase 2 is a poor thermal conductor but good 
electrical conductor. The demand that this sum is maximized sets up a competition between 
the two effective transport properties, and this demand is met by the Schwartz P and 
D structures. We note that for reasons of mathematical analogy, the optimality of these 
bicontinuous composites applies to any of the pair of the following scalar effective properties: 
electrical conductivity, thermal conductivity, dielectric constant, magnetic permeability, and 
diffusion coefficient. 

Since triply periodic minimal surfaces arise in certain cell membranes, it stands to reason 
that such structures must also be relatively stiff elastically in order to withstand mechanical 
loads. The topological property of bicontinuity of these structures suggests that they would 
be mechanically stiff even if one of the phases is a compliant solid or a liquid, provided 
that the other phase is a relatively stiff material. Motivated by this biological example, 
Torquato and Donev asked what are the two-phase composite structures that maximize 
a weighted sum of the dimensionless effective bulk modulus and effective electrical (or ther- 
mal) conductivity. Thus, they set up a competition between a mechanical property and a 
transport property. They demonstrated that triply periodic two-phase bicontinuous com- 
posites with interfaces that are the Schwartz P and D minimal surfaces are again optimal, 
further establishing the multifunctionality of such two-phase systems. 

Based on the aforementioned considerations, it is reasonable to inquire whether the min- 
imal surfaces are optimal for other effective properties, such as flow characteristics. The 
purpose of this paper is to study Stokes (slow viscous) flow through triply porous media 
whose interfaces are the Schwartz P, the Schwartz D, and the Schoen G minimal surfaces. 
In particular, we are interested in computing the fluid permeability fc, which is the key 
macroscopic property for describing Stokes flow through porous media P, 11, Q|. The 
quantity k is the proportionality constant between the average fluid velocity and applied 
pressure gradient in the porous medium as defined by Darcy's law: 

U = --Vp , (1) 

where U is the average fluid velocity, Vp is the applied pressure gradient and fi is the 
dynamic viscosity. Our expectation is that the Schwartz P porous medium will have the 
highest fluid permeability among the three minimal surfaces because it has the smallest 
specific surface (i.e., interface area per unit volume). Surface area alone of course does not 
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necessarily inversely correlate with fluid permeability . We will show that the Schwartz 
P porous medium does indeed have the highest fluid permeability among the three minimal 
surfaces. A natural question to ask is whether the fluid permeability of the Schwartz P porous 
medium is maximal in any sense. Although we do not answer this challenging question here, 
we also compute fluid permeabilities for other triply periodic geometries, such as flow over 
spheres arranged on the sites of a simple lattice and flow inside two different pore-channel 
models at a porosity (volume fraction of pore space) = 1/2. We flnd that the Schwartz 
P porous medium has the highest fluid permeability among all of the six triply periodic 
porous media considered in this paper. We demonstrate that the fluid permeabilities are 
inversely proportional to the corresponding speciflc surfaces for these structures. This leads 
to the reasonable conjecture that the maximal fluid permeability (scaled by the cell length 
squared) for a triply periodic porous medium with a simply connected pore space at a 
porosity = 1/2 is achieved by the structure that globally minimizes the speciflc surface. 

Due to the numerical difficulties in generating complex meshes of various minimal sur- 
faces, the immersed body finite volume method was adopted in this paper . Momen- 
tum forcing is applied on the body surface or inside the body to satisfy the no-slip boundary 
condition on the immersed boundary and also to satisfy the continuity for the each control 
volume element containing the immersed body. This procedure is based on a finite-volume 
approach on a staggered mesh with a fractional step method. The unit domain is chosen to 
be a cube and is digitized into small cubic subvolumes. 

In Section m we review some procedures to represent triply periodic minimal surfaces. 
In Section ITTTj the immersed-boundary finite volume method for fiow through triply periodic 
porous media is described. In Section HV| we verify the accuracy of our numerical methods by 
computing the fiuid permeability for fiow over spheres arranged on the sites of a simple cubic 
and compare them to corresponding analytical/numerical results of Sangani and Acrivos 
[isj l . Here we also present our evaluations of the fiuid permeabilities for other triply periodic 
porous media, including Schwartz P and D surfaces and the Schoen G surface and two 
different pore-channel models. Concluding remarks and a discussion are given in Section IVl 
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II. TRIPLY PERIODIC MINIMAL SURFACES 



Triply periodic minimal surfaces can be characterized exactly using an Enneper- 
Weierstrass (complex integration) representation only for several cases. For Schwartz P 
and D and Schoen G, analytic surfaces are constructed in terms of elliptic functions by the 



Enneper-Weierstrass representation 



2l|, but in practice this is difficult to use nu- 



merically. Recently, it was shown that triply periodic minimal surfaces can be generated as 
the local minima of the scalar order parameter Landau- Ginzburg functional for describing 
ordering phenomena in microemulsion [2^. This Landau- Ginzburg free-energy functional 
2^ can be utilized to calculate numerically a discretization of a potential iplx) such that 
iIj{x) > for points in phase 1 and tplx) < for points in phase 2. At the phase interface, 
iplx) = 0, which in this case is a Schwartz P,D or Schoen G minimal surface. Moreover, the 
minimal surfaces can be approximated by the periodic nodal surface (PNS) of a sum defined 



in terms of the Fourier series l23 
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26l |. because any periodic surface can be expressed 



as the sum of an infinite number Fourier terms. 

In this paper, the Schwartz P and D and the Schoen G surfaces are produced by both the 
Landau- Ginzburg minimization 2^ and PNS approximations with a few Fourier terms 
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26| on a 129^ unit cell. /^From these potentials, one can readily make discretizations of triply 



periodic bicontinuous two-phase composites in which the two-phase interface is a minimal 
surface. Because phase 1 is topologically and geometrically identical to phase 2, either 
phase can be chosen to be the fiuid phase. Indeed, triply periodic minimal surfaces have the 
special property that the fiuid permeability of phase 1 is equal to the fiuid permeability of 
phase 2. Such phase interchange relations are rare in the case of the fiuid permeability. By 
contrast, phase interchange relations in the case of the effective conductivity and effective 
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and 



elastic moduli of two-phase composites have been known for a long time; see 
references therein. We then used an immersed boundary finite-volume method (described 
below) to calculate numerically the fiuid permeabilities of these composites. 
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III. NUMERICAL SIMULATION FOR STOKES FLOW THROUGH TRIPLY PE- 
RIODIC POROUS MEDIA 



We consider the flow field in an infinite porous medium. The Navier-Stokes equations are 
nondimensionalized with t' := fit/{pL'^) and p' := {p — po)/{fiU/L) where t is the time, fi is 
the dynamic viscosity, p is the density, L is the characteristic length, p is the pressure, po is 
the pressure at a reference, U is the characteristic speed. Dropping the primes leads to the 
dimensionless Navier-Stokes equations 



duj 
dt 



+ Re 



duiUj 



dp 
dxi 



+ 



d dui 

doc j doc j 



duj 
dxi 
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where t is the dimensionless time, the dimensionless Cartesian coordinates (j = 1,2 

or 3), Mj's are the corresponding dimensionless velocity components, p is the dimensionless 
pressure and Re = pUL/p is the Reynolds number 
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FIG. 2: Forcing points in staggered mesh. 



In the limit Re — 0, we obtain the dimensionless unsteady Stokes flow equations, which is 
of central concern in this study. These unsteady Stokes equations recover the steady Stokes 
formulation for infinite time advancement. Due to the numerical difficulties in generating 
complex meshes of various minimal surfaces, the immersed boundary finite volume method 

um forcing, /j, is applied to 



was adopted in this paper 



satisfy the no-slip condition on the immersed boundary. The forcing points are located in 
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a staggered fashion like the velocity components defined on a staggered grid (see Fig. 

Therefore, the governing equations for unsteady incompressible stokes flow becomes 

duj _ _dp_ _^ _d_dv^ 
dt dxi dxj dx-i 

^ = 0, 
dxi 

where /j's are the dimensionless momentum forcing components deflned at the cell faces on 
the immersed boundary and inside the body. 

The Stoke s eq uations represented by Eq. (jS)) are advanced with a fully implicit, fractional 
step method [igI . IsoI with the Crank-Nicholson method in convection term. The fractional 
step, or time-split method, approximates the evolution equation by decomposition of the 
operators in it. A four-step time advancement js^ is modifled for time integration of Eq. 
(jni) as follows: 

*" + i^^(«.+<)+/r. (4) 



At dxi 2 dxj dxj 

Ui — iti dp"- 

At dxi ' 

d dp^^^ 1 dui 



dxi dxi At dxi 



(5) 
(6) 
(7) 



At dxi 

Here, the solution -u" at time step n is advanced to "u"'*'^ at time step n + 1 through inter- 
mediate velocities Ui and Ui. Explicit direct forcing scheme for the momentum equation is 
utilized to obtain the momentum forcing value in Eq. (@)) which should be known a priori 



for time advancement 



32] . Discretizing Eq. explicitly in time and rearranging it for fl 



at a forcing point results in 

' At dxi dxj dxj 

where Ui is the velocity imposed at forcing point to enforce = Ui on the immersed 
boundary. To satisfy the no-slip boundary condition on the immersed boundary, Ui is zero 
when the forcing point is at the immersed boundary or inside the body. All the spatial 
derivatives are resolved with the second-order central- difference scheme using a staggered 
mesh system. The discretized momentum equations are factorized into three operators which 
requires inversion of tridiagonal matrices. Time advancement is achieved by sweeping xi-, 



X2-, and X3-direction tridiagonal matrix algorithm at time intervals of n + i, n + |, and 
n + l. 

To drive flow in a unit cell of a triply periodic porous medium, a pressure drop is applied 
from the inlet to the outlet of the channel in one of the three directions (xi in this paper). In 
the other two directions {x2 and X3), standard treatment for velocity and pressure periodic 
boundary conditions are apphed. To solve the elliptic pressure Eq. ©, the Alternating 
Direction Implicit (ADI) method is applied with a relaxation. Non-uniqueness of a pressure 
solution due to applying a pressure drop condition can be resolved with a point constraint 
which imposes an arbitrary pressure value at an arbitrary position. In fact, the pressure 
gradient only is necessary to solve for velocity field. When the triply periodic surface pos- 
sesses simple cubic symmetry, the periodic boundary conditions can be exploited further. 
Because of the mirror image of the flow and pressure field across X2 = 0.5 and X3 = 0.5 in 
the unit size channel simulation, the simple cubic symmetry motivates setting the velocity 
and pressure at the reflection position across the middle cross-section. 

In general, the forcing point of an arbitrary immersed body does not coincide with the 
immersed boundary, but rather inside the body in the non-body fitted Cartesian mesh. This 
requires a velocity interpolation to impose the correct f/j at a forcing point. However, no 
velocity interpolation scheme is necessary to satisfy the no-slip boundary condition on the 
immersed boundary in this paper because in the digitized triply periodic surfaces the forcing 
point of an immersed body always coincides with the immersed boundary. 

IV. NUMERICAL RESULTS 

In this section, we will apply the immersed boundary method to determine Stokes flow 
velocity fields in triply periodic porous media with cubic symmetry. We begin by computing 
the flow field for a periodic array of spheres in order to validate the procedure and the 
numerical results. Subsequently, we apply the numerical procedure to two different pore- 
channel models as well as the bicontinuous structures with interfaces that are the Schwartz P 
and D minimal surfaces and the Schoen G minimal surfaces. The permeability results for all 
six cases are compared at the porosity (p = 1/2- The fluid permeabilities are obtained from 
the imposed pressure drop Ap and average velocity U of the flow field in the cubic periodic 
cell of side length L. Because of the cubic symmetry, the second-order fluid permeability 
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tensor is isotropic (i.e., it is given by kl, where I is the identity tensor) and therefore the 
determination of the permeabihty along the xi axes is sufficient to characterize k. From Eq. 

the fluid permeability k can be made dimensionless by dividing it by L^. Henceforth, 
all permeabilities will be reported in this dimensionless form. 

The computational domain is chosen to be a unit cube and staggered Cartesian meshes 
were generated, depending on the required resolution. All initial velocities and pressures are 
set equal to zero. 



A. Periodic array of spheres 

To verify the numerical accuracy of the immersed boundary method, we compare our 
calculations to the analytical/numerical results of Sangani and Acrivos for flow past 
periodic arrays of spheres centered on the sites of a simple cubic lattice. The immersed 
body (sphere in this case) is constructed using the following criterion. The immersed sphere 
is identifled if the center of each voxel constituting a single sphere centered at (0.5, 0.5, 0.5) 
satisfies the equation 

(xi - 0.5)' + (X2 - 0.5)2 + (a;3 - 0.5)2 < 

where r is the radius of the sphere. As explained in Section ITTTl a pressure drop is applied in 
Xi direction, and standard treatment for velocity periodic boundary conditions are applied 
in X2 and x^ direction. The fluid permeability k obtained from numerical simulations are 
listed in Table |l] for various values of the scaled particle volume fraction x = ic/cmaxY^^, 
where c is the particle volume fraction and Cmax = 7r/6 is the maximum particle volume 
fraction corresponding to the case when the spheres are in contact. Our fluid permeability 
calculations are in good agreement with the results of Sangani and Acrivos |18| . 

Non-uniform (stretching) staggered meshing is used toward the body over the computa- 
tional domain to get the required spatial accuracy for the immersed boundary method, which 
is important for comparison of our permeability data to corresponding analytical/numerical 
results over a smooth sphere (see Fig. EJ. The discretization of a continuous body with a 
flnite resolution obviously affects the numerical accuracy of the computed fluid permeability 
k, which depends on the value of the grid mesh spacing h or the mesh resolution = 1/h. 
We studied the resolution effects on the numerical accuracy at flxed porosity (p = 1/2. The 
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FIG. 3: (Color online) Mesh for flow field over a sphere. 



fluid permeability data and a corresponding linear least-squares fit are shown against 
in Fig. m 

Various mesh resolutions (A^ = 20,40,60,80,128,150,180,200) are used to reduce the 
digitization error over solid regions in these simulations. Generally, fluid permeabilities 
determined from the immersed boundary method for digitized spheres are overestimated 
but this error, as Table HI shows, decreases as the resolution increases. 

B. Bicontinuous Triply Periodic Porous Media 

Before describing our results for the minimal surfaces, we first consider two different pore- 
channel models that are also bicontinuous and triply periodic. The latter is motivated by 
the fact that one-dimensional flow inside channels are characterized by high permeabilities. 
However, since we are interested in macroscopically isotropic permeabilities (i.e., cases when 
the fluid permeability tensor is isotropic), we consider triply periodic channels that intersect 
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FIG. 4: (Color online) Fluid permeability k for Stokes flow past a simple cubic array of equi-sized 
spherical obstacles as a function of for (f> = 1/2. 

in three orthogonal directions (xi, X2 and xs). A pore in this model is the region where the 
channels meet and its geometry depends on the geometry of the channel. In one model, we 
consider square channels and cubic pores. In the second model, we consider circular channels 
and spherical pores (see Fig. EI)- Cross-sections of the two models are shown in Fig. |B1 where 
a and b are the parameters that enable us to control the relative channel/pore geometries. 
For example, in the case of the cubic pore-square channel or spherical pore-circular channel, 
the pore is defined to be minimal when 6 = 0. 

1. Pore- Channel Models 

We calculated the permeability k for three different combinations of a and b at three 
porosities </> = 0.5,0.7 and 0.9 on fine mesh of 128 x 128 x 128. Table ITTl summarizes 
our permeability results for the cubic pore-square channel geometry. Not surprisingly, the 
permeability increases with increasing porosity 0. The minimal pore case (i.e., 6 = 0) 
provided the highest permeability at any porosity. As shown in Fig. [3 channel fiow away 
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Least square fit 
Numerical results 



0.005 



0.01 



0.015 
1/N 



0.02 



0.025 



0.03 



TABLE I: The fluid permeability k for Stokes flow past simple cubic lattice of spheres. Here 
X = {c/cmaxy^^ , where c is the particle volume fraction and Cmax = 7r/6. 



Y 


k (Present work) 


k (Sangani & Acrivos) 


0.1 




X 10^^ 


9 11 07 X 


10-1 


0.2 


84'^fi 


X 10^^ 


3 821 9 X 


10-1 


0.3 




X 10^^ 


2 0805 X 


10-1 


0.4 


1.2516 


X 10^^ 


1 2327 X 


10-1 


0.5 


7 6379 


X 10^2 


7 4668 X 


1 0-2 


0.6 


4.5803 


X 10^2 


4.4501 X 


10-2 


0.7 


2.6059 


X 10-2 


2.5246 X 


10-2 


0.8 


1.3758 


X 10-2 


1.3197 X 


10-2 


0.85 


9.5964 


X 10-3 


9.1516 X 


10-3 


0.9 


6.4803 


X 10-3 


6.1531 X 


10-3 


0.95 


4.2473 


X 10-3 


4.0031 X 


10-3 


1.0 


2.6727 


X 10-3 


2.5203 X 


10-3 



from the channel intersection (both upstream and downstream) has a Poiseuille-type flow 
profile but, in the vicinity of the channel intersection, the flow is no longer confined and 
the associated maximum speed in the X2-X3 plane is reduced. This is illustrated in Fig. 
where the speed profile for the minimal pore case. Increasing the pore size while keeping the 
porosity fixed, leads to an even lower maximum speed, as depicted in Fig. |H1 Consequently, 
the pore region at the channel intersection has a significant effect on the overall flow field 
and hence fluid permeability of the porous media. 

We also computed fluid permeabilities for the spherical pore-circular channel geometry 
but since the trends are the same as in Table ITU we do not present all of these results here. 
In the minimal pore case (6 = 0), the fluid permeability of the circular channel geometry 
{k = 3.4596 X 10-3) is larger than the corresponding square channel case {k = 3.0743 x 10-3). 
Moreover, the circular channel has a lower specific surface than the square channel at the 
same porosity. This point will be discussed further below. 



12 



FIG. 5: (Color online) Unit cells of two different pore-channel models. Left panel: cubic pore-square 
channel. Right panel: spherical pore-circular channel. 
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FIG. 6: Cross-sections of the two different pore-channel models. 
2. Triply Periodic Minimal Surfaces 

Here we report results for the fluid permeability of triply periodic porous media with 
interfaces that are the Schwartz P and D surfaces and the Schoen G surface. The minimal 
surfaces are generated from the aforementioned Landau-Ginzburg potential minimization 
and each has a fine mesh resolution of 129 x 129 x 129. The potential data iIj{x) are 
numerically digitized such that iIj{x) > for points in phase 1 and tplx) < for points in 
phase 2 or vice versa. 

It is found that the computational difficulty due to the geometrical complexities are very 
well treated by the immersed boundary finite volume method. Representative flow flelds 
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TABLE II: The fluid permeability for the cubic pore-square channel model for several porosites 
and geometric parameters a and b. 






a 


b 


k 




0.2500 


0.0000 


3.0744 X 10"^ 


0.5 


0.1324 


0.2500 


0.5310 X 10-^ 




0.0781 


0.3125 


0.0948 X 10-^ 




0.3125 


0.0000 


9.1168 X 10-^ 


0.7 


0.2656 


0.1406 


6.2495 X IQ-^ 




0.2250 


0.2000 


3.7939 X 10-^ 




0.4000 


0.0000 


26.207 X IQ-^ 


0.9 


0.3750 


0.0750 


23.915 X 10-^ 




0.2625 


0.2125 


8.0231 X 10-^ 



over the minimal surfaces are shown in Fig. El Our computed fluid permeabihties for the 
three minimal surfaces are summarized in Table IIIII along with corresponding results for the 
circular and square channels with minimal pores and the simple cubic array of spherical 
obstacles at the porosity (p = 1/2. It should be noted that the Schwartz P surface has 
the highest fluid permeability among these six triply periodic porous media. Included in 
the table is the specific surface s for each of the structures. The fluid permeabilities are 
seen to be inversely proportional to the specific surfaces. Thus, the Schwartz P surface 
has the largest fluid permeability followed by the circular channel, square channel, spherical 
obstacle, and finally the Schoen G and the Schwartz D surface cases. The permeability of 
the circular channel with a minimal pore is actually close in value to that of the Schwartz P 
surface. Evidently, the sharp corner present in the former results in slightly greater energy 
dissipation than in the smooth Schwartz P surface. 
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FIG. 7: (Color) The flow fleld in a unit cell of the triply periodic cubic pore-square channel geometry 
for the minimal pore case (a = 0.250, 6 = 0) at porosity (p = 1/2. Color scheme: red, yellow, green, 
light blue, and blue represent the spectrum of speeds from highest to lowest, respectively. 

TABLE III: The fluid permeability k and speciflc surface s at (/> = 1/2 for six different triply 
periodic porous media models. The pore-channel results are presented for the minimal pore case 
(6 = 0). SC and BCC stand for simple cubic and body centered cubic lattices, respectively. 



Model 


k 


S 


Symmetry 


Schwartz P surface 


3.4765 X 10" 


-3 


2.3705 


SC 


Circular channel/spherical pore 


3.4596 X 10" 


-3 


2.6399 


SC 


Square channel/cubic pore 


3.0743 X IQ- 


-3 


3.0000 


SC 


Spherical obstacle 


3.0591 X 10" 


-3 


3.0780 


SC 


Schoen G surface 


2.2889 X 10" 


-3 


3.1284 


BCC 


Schwartz D surface 


1.4397 X 10- 


-3 


3.9011 


SC 



V. CONCLUSIONS AND DISCUSSION 

We have determined the Stokes flow fields and associated fiuid permeabilities through six 
difi'erent triply porous media: the Schwartz P, the Schwartz D, and the Schoen G minimal 
surfaces, two different pore-channel models, and a simple cubic array of spherical obstacles. 
The latter model served as a benchmark to verify the accuracy of our numerical procedure: 
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FIG. 8: (Color) The flow fleld in a unit cell of the triply periodic cubic pore-square channel 
geometry for the case a = 0.132, b = 0.250 at porosity (p = 1/2. The color scheme for speeds is the 
same as described in Fig. [71 




FIG. 9: (Color) Flow fields for triply periodic porous media with interfaces that are minimal 
surfaces. Left panel: Schwartz P surface. Center panel: Schwartz D surface. Right panel: Schoen 
G surface. The color scheme for speeds is the same as described in Fig. [7| 

the immersed boundary finite volume method. We found good agreement with previous 
results for the simple cubic array. Interestingly, we found that the Schwartz P porous 
medium has the highest fiuid permeability among all of the six triply periodic porous media 
considered in this paper. The fiuid permeabilities were shown to be inversely proportional 
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to the specific surfaces for these bicontinuous structures. This might lead one to conjecture 
that the maximal fluid permeability (scaled by the cell length squared) for a triply periodic 
porous medium at a porosity = 1/2 is achieved by the structure that globally minimizes 
the specific surface. As we pointed out in the Introduction, the specific surface does not 
necessarily correlate with the permeability. However, the pore spaces of structures with 
large permeabilities are expected to be simply connected and therefore it would not be 
unreasonable for the permeability to be inversely proportional to the specific surface in 
these instances. 

The determination of the structure that globally optimizes the fluid permeability is a 
highly challenging problem. It is only very recently that optimal bounds on k have been 
identified. By contrast, optimal bounds on the effective conductivity and elastic moduli of 
two-phase composite media have long been known 
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In particular, Torquato and 



Pham j33| showed that the so-called void upper bounds on k for flow in anisotropic media 
with a single preferred direction are optimal for certain coated-cylinders models. Such 
structures arise in practice in sea ice, for example, which has a microstructure that is well 
modeled by oriented cylindrical channels Q]. However, these bounds are not applicable for 
our present purposes because we would like to ascertain the triply periodic structure that 
maximizes the isotropic permeability. Even if the conjecture that the maximal dimensionless 
fluid permeability for a triply periodic porous medium with a simply connected pore space at 
a porosity (p = 1/2 is achieved by the structure that globally minimizes the specific surface 
is true, finding the minimizing specific surface structure is highly nontrivial. For example, 
Kelvin's problem, the determination of the space-filling arrangement of congruent closed 
cells of equal volume that minimizes the surface area, is still an open question, although it 



is believed that the Weaire-Phelan structure 



is a good candidate. In a future work, we 



will determine whether the Schwartz P triply periodic porous medium is a locally maximal 
structure for the isotropic fluid permeability. 
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